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O Abstract 



We study fracture processes within a stochastic fiber-bundle model where it is assumed that after the failure of a fiber, each intact 
fiber obtains a random fraction of the failing load. Within a Markov approximation, the breakdown properties of this model can 
!_i be reduced to the solution of an integral equation. As examples we consider two different versions of this model that both can 
, ^ interpolate between global and local load redistribution. For the strength thresholds of the individual fibers, we consider a Weibull 
distribution and a uniform distribution, both truncated below a given initial stress. The breakdown behavior of our models is 
compared with corresponding results of other fiber-bundle models. 
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■j. 1. Introduction 



Fracture processes in heterogeneous materials are an impor- 
tant technological problem that has attracted the interest of 
the scientific community since a long time [1—4]. Due to the 
complex interaction between failures and the subsequent re- 
distribution of local stresses, the development of adequate sta- 
tistical models for fracture propagation is an extremely hard 
and challenging undertaking. The probably most important 
class of approaches to the study of fracture processes is that of 
fiber-bundle models (FBM's) [5-15]. Despite their simplicity, 
FBM's are able to describe the main processes that can lead to 
a propagation of fractures and eventually to a complete break- 
down of real heterogeneous materials. 

Fiber-bundle models refer to a bundle of N parallel fibers that 
are clamped at both ends and stretched by a common force F. 
The fibers have a stochastic distribution of individual strength 
thresholds, and the different versions of FBM's that have been 
considered can be distinguished by their assumptions with re- 
spect to the stress redistribution after the failure of one of the 
fibers. The usual experimental setup considered and analyzed 
in the FBM literature can be described as follows: The force F 
is gradually increased from zero until the weakest fiber breaks, 
and the transfer of its stress to the surviving fibers may then 
induce an avalanche of subsequent failures. If the fiber bundle 
reaches an equilibrium with no further failures, the force F is 
increased again until the next fiber breaks, and this procedure 
is repeated up to the complete breakdown of the entire bun- 
dle. The main quantities of interest in connection with this pro- 
cedure are the distribution of avalanche sizes and the ultimate 
strength of the fiber bundle, defined as the maximum stress F/N 
the system can support before it breaks down completely. An 
alternative but equivalent procedure is to apply a finite force F 
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to the system, so that immediately all fibers with a strength 
threshold smaller than F/N fail. The ultimate strength of the 
fiber bundle is then determined by the maximum F/N value 
that does not lead to a failure of the entire system. 

The oldest and most well-known FBM is that where the 
stress of a failing fiber is distributed equally between the sur- 
viving fibers (global load sharing, GLS) [5]. For this model, the 
strength of the fiber bundle as well as the form of the avalanche- 
size distribution can be determined analytically [6, 7]. 

Local load sharing (LLS) fiber-bundle models, on the other 
hand, are much more difficult to analyze [6, 8, 9]. In these 
models, the stress of a failing fiber is only transferred locally, 
typically to the surviving nearest neighbors. LLS models have 
been studied mainly via Monte-Carlo simulations and analyti- 
cal results have only been obtained for one-dimensional models 
with essentially nearest-neighbor load transfer. 

In most studies, the strength thresholds of the individual 
fibers are assumed to be distributed according to a Weibull dis- 
tribution (typically with a Weibull index k = 2). For simplicity, 
however, uniform strength-threshold distributions (sometimes 
with a finite lower cutoff) have also been considered [6, 7, 11]. 

The two idealized extremes of global load sharing and of load 
transfer to nearest neighbors only are not adequate assumptions 
for most real systems. Attempts have therefore been made to 
interpolate between GLS and LLS behavior [10-13]. Hidalgo 
et al. [12], e.g., assume that the stress transfer after a failure de- 
cays as l/r 7 , where r is the distance from the broken fiber, and 
they study the failure-propagation behavior of such a system as 
a function of y. A similar model (with y — 3) has been studied 
by Curtin [14]. 

In Ref. [16], we have introduced a class of failure- 
propagation models that can represent, in a stochastic sense, 
the main characteristics of realistic load-redistribution mecha- 
nisms, but are still amenable to an analytical treatment. We 
have applied our approach to an illustrative prototype exam- 
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pie for cascading failure propagation in large infrastructure net- 
works, e.g., power grids. In particular, we have analyzed the 
probability of a complete system breakdown after an initial 
failure, and we have found that the model exhibits interesting 
critical dependencies on parameters that characterize the fail- 
ure tolerance of the individual elements and the range of load 
redistribution after a failure. 

In this paper, we apply our stochastic approach to the prob- 
lem of fracture propagation in fiber bundles and analyze two 
models that interpolate between global and local stress redis- 
tribution. The first is a stochastic version of the r~ y -model of 
Hidalgo et al. [12], and the second a fiber-bundle equivalent to 
our prototype example of Ref. [16]. We consider an experimen- 
tal setup where initially all fibers carry the same stress cr and 
where the individual strength thresholds cr tn are randomly dis- 
tributed according to a probability density p t h(o~±) that is zero 
for cr± < ctq. We then examine the consequences of the break- 
ing of a single fiber and concentrate on the calculation of the 
following quantities: 

(i) the no-cascade probability P nc (cr ), i.e., the probability 
that an initial failure does not induce any further failures; 

(ii) the breakdown probability P\,(cr Q ), i.e., the probability 
that an initial failure leads to a breakdown of the entire 
fiber bundle; and 

(iii) the critical stress cr^, defined as the largest cr* such that 
^b(o"o) = for all cr < cr*. 

We note that, in contrast to most of the other fiber-bundle 
models, our stochastic models neglect any spatial correlations 
in the load transfer after a failure. The only other model, as 
far as we are aware of, that also uses a stochastic (rather than 
a spatially correlated) stress redistribution is that of Dalton et 
al. [15], where it is assumed that the load of a failing fiber is 
transferred to a fixed number (n = 1, 2, . . . ) of randomly chosen 
surviving fibers. 

In addition, we use strength-threshold distributions that are 
truncated below the initially applied stress cr — in contrast to 
most models studied in the fiber-bundle literature. There exist, 
however, a number of investigations that also consider truncated 
threshold distributions [7, 11, 4], so that a direct comparison 
with our results can be made. 

In Sect. 2, we introduce our stochastic load-redistribution 
model and describe its application to fracture processes. Sub- 
sequently, in Sect. 3, we describe a Markov approximation of 
the model which leads to a description in terms of generalized 
branching processes. In Sects. 4 and 5, we analyze the two dif- 
ferent model variants mentioned above, and final conclusions 
are given in Sect. 6. 

2. Stochastic load-redistribution model 

2.1. Load-redistribution rule and cascade model 

We shall consider a bundle consisting of N fibers subjected 
to an external force F. We assume that the initial stress 
cr = F/N of all fibers is equal and that the strength thresh- 
olds cr th of the individual fibers are randomly distributed ac- 
cording to a probability density /?th(o"th)- For our setup of 



fracture-propagation experiments, we assume that we start from 
a finite stress cr > and that all fibers with thresholds smaller 
than cr have been removed [4]. Thus, the threshold distribution 
has to fulfill /?th(o"th) = for all cr th < cr . On the other hand, 
we are interested in a situation where already an infinitesimal 
increase of the external force leads — with probability one in the 
limit N — > oo — to the breaking of exactly one fiber, and we thus 
require p t h(o-±=o- ) > 0. 

When a fiber breaks, its stress has to be taken over by the 
remaining intact fibers of the bundle. In our stochastic load- 
redistribution model [16], we assume that this process can be 
described by a rule of the form 

cr -> cr' = cr + cr f A . (1) 

Here, cr (cr') is the stress of an intact fiber before (after) the 
failure of a fiber with stress o~f , and the load-redistribution fac- 
tor A is a random number drawn independently from the same 
distribution /?a(A) for each of the intact fibers. Note that for the 
initial failure, both cr and cr f are given by the initial stress cr . In 
the special case A = l/(N - 1) of a uniform load-redistribution, 
the form (1) reduces exactly to a GLS rule. For a general non- 
uniform load redistribution, we require that the failed stress will 
be shared on average by the remaining intact elements, i.e., the 
mean of A has to fulfill the condition 



Due to the stress increment a fiber has obtained after a failure 
event, its stress itself might be above its critical threshold. In 
general, the initially failing fiber might thus induce the failure 
of N f ' > 1 other fibers, thereby starting a failure cascade. We 
then assume that all overloaded fibers fail simultaneously and 
that their stress is again redistributed according to the rule (1), 
where now both the pre-failure stresses cr of each of the N-N^ l) 
intact fibers and the stresses cr f of each of the N ( f r> failing fibers 
will, in general, be different random variables. If this process 

(2) 

leads to the further overloading of N f > 1 fibers, it continues 
to a next cascade stage, and so on. Eventually, either the bundle 
stabilizes again, i.e., all fibers are stressed below their respec- 
tive strength thresholds, or it breaks down completely, i.e., all 
N fibers fail. 

Note that, in general, during the failure cascade, the load 
redistribution and hence the A-distribution will be modified. 
Whereas the details of such a modification, which can depend 
on topological changes, are very difficult to model, one at least 
has to take into account one dominant effect: As the number 
of fibers that are still intact at cascade stage s decreases, 

in ° ' 

the mean (A) has to increase in accordance with Eq. (2) with 
N replaced by A^ n ?) . Below, we will discuss how to fulfill this 
requirement for the chosen forms of load redistribution. 

2.2. Truncated strength-threshold distributions 

For our setup of fracture-propagation experiments, we have 
to truncate the distribution of strength thresholds below cr^ = 
cr . In the literature on fiber-bundle models, the strength thresh- 
olds crth of the individual fibers are usually assumed to be 
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distributed according to a Weibull distribution with density 
k cr^' 1 exp(-cr^), where mostly the Weibull index k = 2 is used. 
Truncation then leads to a distribution of the form 



Pth(o-< b ) = r lh 



k o-f 1 exp(cr* - erf ) if cr < cr th 



else. 



(3) 



In addition, we shall also consider uniform distributions that are 
truncated below the initial stress cr : 

1 

— it iTn < (T,u < 1 

(4) 



Pth(o-th) = <l-(To 





if cr < cr th < 1 

else. 



In Sects. 4 and 5, we summarize and discuss the correspond- 
ing results for P nc (o"o), Pb(o"o) an d °"o ^ or two different load- 
redistribution models and for the two threshold distributions (3) 
and (4). 

3. Generalized-branching-process approximation 

The dynamics of the stochastic cascade model described in 
the previous section and the quantities P n c(o"o)> Pb(o~o) and cr^ 
can only be obtained exactly by means of Monte-Carlo simula- 
tions. In the limit of large system sizes N — > oo, however, we 
can achieve an approximate description of the cascade dynam- 
ics by noting the following points: 

(i) The failure of a fiber leaves the stress in the majority of 
the intact fibers nearly unaffected, i.e., maximally leads 
to changes of the order of 1 /N. Thus, along the failure 
cascade, the stress of the intact fibers is approximately 
given by the initial stress cr . 

(ii) The remaining number of fibers always stays infinitely 
large and thus the number of further failures induced by a 
failing fiber is distributed according to a Poisson distribu- 
tion. 

(iii) The interaction between different failures can be ne- 
glected, i.e., in the case of several induced failures, the 
failure cascades resulting from each of these failures can 
be treated as being independent. 

Under these assumptions, the cascade dynamics becomes 
Markovian [18] if we choose the point process of the stresses 
of the failed fibers as underlying state space. This point pro- 
cess on the semi-infinite interval [ctq, oo) is independent [17], 
and the failure dynamics can thus be described by a generalized 
branching process [19] with characteristic functional 

G[u; cr { ] = exp {ju(^f) J" > a-*; cr f ) - 1 } 

(5) 

for the point process induced by a single failure with stress <x f . 
Here, p(o~' f \(r' f > cr^; cr f ) denotes the conditional probability 
density that the induced failure resulting from the breaking of a 
fiber with stress cr f occurs with a stress crj.. For given distribu- 
tions of the load-redistribution factors A and the critical thresh- 
olds cr±, this quantity can be readily calculated from Eq. (1). 
This also holds true for the mean number of failures, 



induced by the breaking of a fiber with stress <x f . We remark 
that in order for this quantity to be finite in the limit N -» oo, 
the probability P(cr' { > cr±\crf) for the induced failure of a given 
intact fiber has to vanish as 1 IN. This is in accordance with 
the requirement (2) for the mean of the load-redistribution fac- 
tors A. Finally, in Eq. (5), u denotes an arbitrary non-negative 
test function on the interval [cr , oo). 

For later use, we note that from the mean number of induced 
failures, one directly obtains the no-cascade probability, which 
is given by 



^nc(o-o) = [1 - P(cr' f > o-thlo-o)] 
In the limit N — > oo, this relation becomes 

fnc(o-o) = exp [ - yU(<T )] . 



N-l 



(7) 



(8) 



In principle, the properties of the later cascade stages and 
thus the full cascade dynamics can be obtained in a recursive 
way from the functional (5) (see Ref. [19]). Here, we are only 
interested in the question whether an initial failure leads to the 
breakdown of the entire fiber bundle. It can be shown [19] that 
this question can be answered by solving the integral equation 

1 - Pb(o-f) = exp | - n{(Tf) J dcr' t p((r' f \cr' { > cr th ; cr f ) 

x^K)} (9) 

for the probability Pb(o"f) that an initial failure with stress <Xf 
leads to the breakdown of the entire bundle. Using Eq. (6) 
and the load-redistribution rule (1), together with the fact that, 
within the approximation considered, the pre-failure stress <x is 
just given by the initial stress cr , we can rewrite the integral 
equation (9) in the form 

CO 

1 - P b (cr f ) = exp { - i- J d^ J F th (^)/5 A (^^) 



xPbK)}, 



(10) 



where /5a(A) = limj V _, 00 Np&(A). Here, we have also used the 
relation 



p(o-' f \o-' f > cr th ; cr f ) P(o-' { > cr th |cr f ) = p(cr' f ; cr' f > cr th |cr f ) (11) 



and introduced the cumulative distribution function Fth(c t h) = 
dcr' lh p th (cr' lh ) corresponding to the threshold distribution 
Pth(o"th). Under quite general assumptions [19], the unique 
solution of the integral equations (9) or (10) can be found by 
means of an iterative procedure starting from an arbitrary ini- 
tial guess for Pb(cr f ). The so obtained probability function can 
then be evaluated at the initial stress <tq to obtain the proba- 
bility Pb(o~o) for the breakdown of a fiber bundle in the setup 
described in Sect. 1 . 

We finally remark that the interpretation of Eq. (9) becomes 
clear if one writes the exponential on the right-hand side in the 
form 



Mo- f ) = (AT-l)P( ( r' f XTfrlo-f), 



(6) 



1 



-Pb(o-f) = lim 1- do-' { p(o-' { ;o-' { >o- th \o-f)P b (o-' f ) 



(12) 



where we have again used Eqs. (6) and (1 1). Hence, the proba- 
bility that no breakdown occurs after the failure of a fiber with 
stress cr f is equal — in the limit N — > oo — to the probability that 
none of the induced failures with stress cr' f leads to a breakdown. 



4. Stochastic model for range-dependent load redistribu- 
tion (y-model) 



Hidalgo et al. [12] proposed a FBM where the stress trans- 
fer after a failure decays with the distance r between the failing 
fiber and the one affected by the failure as a power law Z/r r . 
Here, Z is a normalization factor which ensures that the total 
load is conserved. They furthermore assumed that all fibers are 
arranged on a two-dimensional square lattice. By varying the 
exponent y > 0, they were then able to study the transition be- 
tween a GLS rule (for y — > 0) and a LLS rule (for y — > oo). Note 
that due to the infinite range of the power-law transfer function, 
the latter situation of a strictly local load-transfer can only be 
achieved in an approximate sense. For the case of Weibull- 
distributed strength thresholds, a Monte-Carlo analysis of this 
range-dependent load-transfer model showed that for an expo- 
nent y < 2, the model behaves essentially as a FBM with GLS 
rule and, in particular, a finite critical stress value was observed. 
For larger y, there is a transition to the LLS case with a critical 
stress that vanishes in the large system-size limit. In Ref. [11], 
the same model has been analyzed for uniform threshold dis- 
tributions with a lower cutoff cr L , and in this case, the critical 
stress remains finite for all values of y (if cr L > 0). 

In the following, we study a stochastic version of this model, 
which we shall call "y-model". It is based on the assumption 
that the position of the fibers is uniformly distributed within 
the two-dimensional cross-section of the bundle. Upon failure 
of a fiber we then randomly pick the affected fibers from this 
uniform distribution and calculate the load-transfer factor A ac- 
cording to the (random) distance. We will now first derive the 
corresponding A-distribution, then analyze the properties of the 
resulting model, and finally compare its results to the ones ob- 
tained in Ref. [11]. 



4.1. Distribution of load-redistribution factors 

For reasons of simplicity, we assume that the broken fiber 
is in the center of a hollow cylinder with inner (outer) ra- 
dius r m ; n (r max ) containing N - \ intact fibers uniformly dis- 
tributed with area density g in the cross-sectional area of size 



A = n(r n 



'"min )■ Note that in contrast to Refs. [11, 12], we 



consider a uniform distribution of the fiber positions and thus 
have to introduce a lower cut-off for the distance r to prevent a 
divergence at small distances. 

For this uniform spatial distribution and the given distance 
dependence of the load transfer, we can then readily derive the 
probability distribution for the load-redistribution factors A ap- 



pearing in Eq. (1): 

'"max 

p A (A):= ^ J drrdiA-Z/r?) 



(13) 



2 r+2 

2nZy A r 



Ay 







for A min < A < A max 
otherwise. 



(14) 



Here, the lower and upper cutoffs for A are given by A„ 



Z/r max y and A n 



Z/r mm y , respectively. In the stochastic 



model, we fix the constant Z by imposing Eq. (2), which states 
that, on average, the load is redistributed to the remaining ele- 
ments. This yields the normalization constant 



2 _ 2 
f max — f min 



2-y 

2(N - 1) r max 2 -r - r n 



(15) 



for y + 2; the special case y = 2 can be readily treated by 
considering the limit y — » 2. 

We can write Eqs. (14) and (15) in a more convenient form 
by introducing a dimensionless length L = r max /r mm . This first 
allows us to express the number of intact fibers as 



N- 1 = s(L 2 - l) , 



(16) 



where s = nr m i n 2 g « 1 is the average number of fibers in the 
vicinity of the failing fiber. From Eq. (16), we obtain, e.g., 
s = it I A if we assume that the model describes the continuous 
approximation of fibers located on a quadratic lattice (with lat- 
tice constant r min ) consisting ofN~ n(L/2) 2 sites inside a circle 
of the given radius. Furthermore, we can write the A-cutoffs in 
the form 



A m j n — _ _„ 



11 2-y 



Amax — _ 



2s L 2 1 - V- 1 
1 2-y 



2s Lr- 2 - 1 



(17) 



(18) 



where again the case y = 2 has to be treated as limit. The 
probability distribution is then given by the power-law form 



2 1 1 

PA(A) = y Z^T A~ 

Y L 1 "ma 



7+2 

(aH 7 (19) 

\ "max / 



for A min < A < A max . 

For later use, we note that in the limit of large system sizes, 
the lower A-cutoff always scales to zero; A mm — > for L — > oo. 
The behavior of the upper cutoff, however, strongly depends on 
the exponent y: For y < 2, A max also vanishes in the limit L — > 
oo, whereas for y > 2, it converges to a finite value: A max — > 
(y-2)/2sforL-^ oo. 

4.2. Mean number of directly induced fiber failures 

We now calculate the mean number //(<Xf) of fiber failures 
resulting directly from the breaking of a fiber with stress cr f . 
From rule (1) and Eq. (6), we obtain 



/i(tr f ) = (N - l)P(cr + o-f A > o- m ) 



(JV-1) dA/7 A (A)F m (o-o + cr f A). 



(20) 



With Eqs. (16)— (19), we can write this expression in the form 

1 



— f dxx 

y J 



y+2 



°~f ^max 

x) . (21) 



For the evaluation of this integral in the limit L — > oo, where 
the integrand becomes singular, it is useful to consider the cases 
y < 2 and y > 2 separately. 

(i) y < 2: As mentioned above, we then have A max -> Ofor 
L — > oo and thus consider the Taylor expansion 



Fth(o- () + cr f A max x) = ^ 



(o-f A max x) 1 



F«Vo), (22) 



/=i 



where we have used that F t h(o"o) = 0. Inserting this expansion 
into Eq. (21), we find with Eq. (18) that only the first order 
term 1=1 leads to a non-vanishing contribution in the limit L — > 
oo. This yields 

fi(o-f) -> cr f pfrio-o) . (23) 

Note that this result is independent of the exponent y. 

(ii) y > 2: In this case, we can insert the finite asymptotic 
value for A max , A max = (y - 2)/2s for L — > oo, into Eq. (21) and 
write the mean number of failures as the improper integral 

l 

2s r _zi2 

yu(cr f ) = — dxx y F th (o- + o-f(y-2)/(2s)x). (24) 
o 

Here, the convergence of the integral at the lower boundary is 
guaranteed since, because the strength-threshold distribution is 
truncated below cr , we have /*\h(o"o + o"f (y-2)/(2s)x) = O(x). 

4.3. Breakdown probability 

For the evaluation of the breakdown probability, we have to 
solve the integral equation (10), which in the present case as- 
sumes the form 



1 - Pb(o-f) = exp { - 



2s 



o"o+o"f A m 



yo-f A n 



- cr 



xFthK)P b K)J. 

Alternatively, we can write this equation as 

l 



(25) 



1 - /y<r ( ) = expj - y f dx x ' /7 >ii , <r " : <^A max .v) 

X P b (cr + x) . (26) 



The integral on the right-hand side of this equation has to be 
evaluated in the limit of infinitely large systems (L — > oo). 
Again, we treat the cases y < 2 and y > 2 separately. 

(i) y < 2: As above, in Eq. (22), we expand the distribu- 
tion function f\h(o"th) around cr th = cr and now furthermore 



assume that such an expansion is also valid for the breakdown 
probability, 



P b (o- Q + o- f A max x) = ^ 



(c f A max xf 



m=0 



ml 



ffVo). 



(27) 



Inserting these expansions into the integral in Eq. (26), we find 
with Eq. (18) that the various terms behave as £ 2(1 ~'~ m) for large 
L. Thus, only the lowest order terms / = 1 and m = survive 
in this limit. The integral equation (26) hence simplifies to the 
transcendental equation 



1 - Pb(o-f) = exp [ - o-f pth(cro) Pb(o-o)] , 



(28) 



which, again, is y-independent. In this regime, the critical 
stress cTq is given by the condition that the mean number of 
directly induced failures equals unity: 



McT()) = crSAhK))= 1 



(29) 



(ii) y > 2: As for the evaluation of the mean number of 
induced failures, cf. Eq. (24), we use the asymptotic value of 
A max and replace the integral in the limit L — > oo by an improper 
one [cf. also remark after Eq. (24)]. This leads to the integral 
equation 

l 

1 - Pb(o- f ) = exp | - y J dx x^-TFUo-o + cr f (y - 2)/(2s)x) 
o 



x P b (cr () + o- f (y - 2)/(2s)x)\ , (30) 

which, in general, can only be solved numerically, e.g., by 
means of an iterative procedure. In particular, the critical 
stress o"q is not determined by a simple relation like Eq. (29) 
but has to be determined from the full solution of Eq. (30). 

4.4. Results 

4.4.1. Uniform distribution of strength thresholds 

In the case of a uniform distribution (4) of the strength thresh- 
olds, the mean number of induced failures (21) can be evaluated 
explicitly. With cr f = cr we then obtain from Eq. (8) the no- 
cascade probability 



Pnc(o-o) = exp 



\ l-crj 



if y < 2 or o-q < [1 + (y - 2)/(2s)]-\ and 



Pnc(o"o) = exp 



-I 



y-2 

2s 



n 



1 - (TO 



2/7 



(31a) 



(31b) 



otherwise. For the calculation of the breakdown probabil- 
ity Pb(o"o), the transcendental equation (28) (for y < 2) or the 
integral equation (30) (for y > 2) have to be solved numerically. 

In Fig. 1 we show the no-cascade probability P n c(o"o) and 
the breakdown probability Pb(o"o) as a function of the initial 
stress o-q for different values of the exponent y in Eq. (19). The 
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Figure 1 : y-model with s = n/4 and uniform distribution of strength thresholds, 
(a) No-cascade probability P nc and (b) breakdown probability Pb as a function 
of the initial stress <tq for different values of the exponents y. Lines: Results 
from Eq. (31) [panel (a)] and Eqs. (28) and (30) [panel (b)], respectively. Sym- 
bols: Results from Monte-Carlo simulations of the failure process for L = 64 
averaged over 1000 realizations. The statistical error is of the order of the size 
of the symbols. 



approximate results within the Markov approximation (lines) 
are compared with Monte-Carlo simulations (symbols) of the 
failure dynamics generated by the load redistribution (1). We 
note that in the Monte-Carlo simulations, we have to ensure 
that condition (2) stays fulfilled during the entire cascade pro- 
cess. This is done by replacing L in Eqs. (17)— (19) by L e s = 
y/Ni n /s + 1, where denotes the number of remaining intact 
fibers. 

With increasing initial stress <x , we observe a gradual de- 
crease of the no-cascade probability from one to zero [cf. 
Fig. 1(a)]. In contrast, the breakdown probability [Fig. 1(b)] 
exhibits a critical behavior: There is a y-dependent critical 
stress crS such that for <x < crS, the probability of a break- 
down of the fiber bundle vanishes exactly. The dependence 
of the critical stress on the exponent y is illustrated in Fig. 2. 
The value <x£ = 1/2 for y < 2 can be readily determined from 
Eqs. (4) and (29) and exactly reproduces the value of the FBM 
with GLS rule [10]. For y > 2, and hence smaller effective 
"range" of the stress redistribution, we first observe a transition 
to a regime, where the critical stress decreases with increasing y 
down to a minimal value crjc ~ 0.26. For even larger y, y > 7, 
the critical stress increases again. We remark, however, that for 
large y, it becomes numerically rather difficult to find the pre- 
cise location of the critical transition because the onset of the 
regime with a finite breakdown probability becomes more and 
more flat [cf. Fig. 1(b)]. 
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Figure 2: y-model with j- = jr/4. Critical stress d"q as a function of the ex- 
ponent y. Solid line: uniformly distributed strength thresholds. Dashed line: 
Weibull-distributed (k = 2) strength thresholds. 



It is interesting to compare our results with the ones obtained 
for the variable -range load-redistribution model of Ref. [1 1], in 
particular, the behavior for the case of the failure stress being 
equal to the cutoff strength [cf. Fig. 1(a) of Ref. [11]]. In both 
models, we observe a critical value of y sa 2 above which a tran- 
sition from a GLS regime to one with short-ranged stress trans- 
fer and smaller cr^-value takes place. Within our model, we can 
trace back this transition to a change in the load-redistribution 
distribution, in particular, the asymptotic value of A max . Fur- 
thermore, the critical stress values of both models agree rather 
well up to y < 7. For even larger exponents y, we find an 
increase of the critical stress, which cannot be observed in the 
more microscopic model of Ref. [11]. This discrepancy proba- 
bly results from a breakdown of the continuum approximation 
upon which the distribution (13) is based. The deficiency of 
this approximation for large y is also reflected in the fact that 
for y > 2 + 2s, the asymptotic value for A max becomes larger 
than unity, which means that a single fiber may receive a stress 
increment that is higher than the stress of the failing fiber. In 
order to prevent such a pathological behavior, a more sophisti- 
cated load-redistribution model has to be used. 

We finally note that the results from the Markov approxima- 
tion, i.e., the generalized branching process description, agree 
very well with the ones obtained by Monte-Carlo simulations 
of the failure process. Around the critical transition, some de- 
viations for the breakdown probability can be observed, which, 
however, decrease with increasing system size N and, thus, rep- 
resent finite-size effects [16]. 

4.4.2. Weibull distribution of strength thresholds 

For the case of the truncated Weibull distribution (3) of 
strength thresholds, we have to evaluate for y > 2 both the no- 
cascade probability and the breakdown probability numerically 
from Eqs. (24) and (30), respectively. 

The results as a function of the initial stress ctq are depicted 
in Fig. 3, where we have chosen here and in the following a 
Weibull index of k — 2. Comparing with the case of a uniform 
distribution of the strength thresholds, we find qualitatively the 
same behavior. In particular, we identify a critical transition at 
a y-dependent stress <x£ (cf. Fig. 2), and again, for y < 2, the 
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Figure 3: y-model with .v = ;r/4 and Weibull-distributed strength thresholds 
(k = 2). (a) No-cascade probability P nc and (b) breakdown probability Pb as a 
function of the initial stress <tq for different values of the exponents y. Lines: 
Results from Eqs. (23) and (24) [panel (a)] and Eqs. (28) and (30) [panel (b)], 
respectively. Symbols: Results from Monte-Carlo simulations of the failure 
process for L = 64 averaged over 1000 realizations. The statistical error is of 
the order of the size of the symbols. 



result cTq = for the GLS case [5] is recovered exactly 

from Eq. (29) together with the distribution (3). 

Figure 4 shows the dependence of the breakdown probabil- 
ity on the exponent y for fixed initial stress cr () . In accordance 
with Eq. (28), the breakdown probability is y-independent for 
y < 2 and assumes the GLS value. In the case y > 2, we 
find for <x > 1 a regime with a monotonic decrease of fb to- 
wards zero as a function of y. For smaller ctq, the breakdown 
probability assumes a maximum at a certain y-value and then 
decreases again towards zero. Finally, for <tq smaller than the 
critical stress of the GLS model but larger than the minimal 
stress observed in Fig. 2, an increase of y eventually leads to a 
destabilization of the system, i.e., a non-vanishing breakdown 
probability, above some critical y-value. 

5. Simple bimodal load-redistribution model (Ao-model) 

In Ref. [16], we have introduced a simple prototype model 
that interpolates between the limiting cases of global load re- 
distribution and the transfer of the failing load to a single other 
element. 

The model, which we shall call "Ao-model", is characterized 
by a bimodal distribution of the load-redistribution factors A, 



A = 



I A(> with probability po 

1 with probability 1 - po , 



(32) 
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Figure 4: y-model with s = tt/4 and Weibull-distributed strength thresholds 
(k = 2). Breakdown probability Pb as a function of the exponent y for different 
values of the initial stress <tq. Lines: Results from Eqs. (28) and (30). Sym- 
bols: Results from Monte-Carlo simulations of the failure process for L = 64 
averaged over 1000 realizations. The statistical error is of the order of the size 
of the symbols. 



i.e., after the failure of an element with stress cry, the stress <x 
of a still intact element is increased to <x' = cr + Aq <x f with 
probability po and remains unchanged with probability 1 - po. 

We further require that the sum of the induced stress incre- 
ments is, on average, equal to the stress of the failing element 
and that A < 1 . It follows that 



Po A 



1 



1 



N — 1 N-l 



< Ao < 1 



(33) 



Ao = l/(N - 1) then corresponds to the limiting case of global 
stress redistribution and Ao = 1 to the case where the failing 
load is transferred, on average, to a single other element. 

The probability that after the failure of a fiber with stress <Xf, 
a still intact fiber also fails can be written as 



Pq P(cr th < cr + A cr f ) = 



1 



(N- 1)A 

and the mean number of induced failures becomes 
1 



P(o~ th < o"o + Ao<Xf ) (34) 



M(o"f) 



P(cr th < cr + A Of) . 



(35) 



In these expressions, we have neglected that a fraction of the 
still intact fibers at later cascade stages may carry a stress larger 
than o"o. It can be shown, however, that this is a finite-size 
effect, i.e., this fraction vanishes as — > oo [16]. 

The no-cascade probability then follows directly from Eq. (8) 
by using Eq. (35) with <Tf = ctq: 



fnc(o-o) = exp 



1 



•P(cr th <cr (l + A )) 



(36) 



To calculate the breakdown probability Pb(o"oX we use Eq. (10), 
which reduces in the present case to the recursion relation 



P b (o-„) = 1 - exp[-ju(cr„)P b (cr„ + i)] , 



where 



cr„ = <x - 



1 - A" +1 
1 - Ao 



n = 0, 1, 



(37) 



(38) 
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This recursion can be solved numerically to an arbitrary degree 
of accuracy by starting at a high enough value of n, say n s , and 
setting Pb(n s ) = 1 . 

Finally, it can be shown that the critical stress <fL is deter- 
mined by 

O-Q 



fiio-n^cc) = 1 , CT„- 



1 - An 



i.e., by 



5.7. Results 



1 / A cr \ , 



(39) 



(40) 



5.1.1. Uniform distribution of strength thresholds 
For a uniform threshold distribution (4) we obtain 

/ AnCTf \ 

P(o- c < cr + A o- f ) = mm 1, , 

V 1 - cr n / 



(41) 



and it follows that 



Pnc(o-O) 



exp ( _ i) 



otherwise 



and 



1- Ap 

2- A 



(42) 



(43) 



The breakdown probability Pb(o~o) is determined by numeri- 
cally solving the recursion defined in Eqs. (37) and (38) with 

Mcr^minfl,^). (44) 
A \ 1 - cr / 

5.7.2. Weibull distribution of strength thresholds 

For the case of a truncated Weibull distribution of strength 
thresholds, Eq. (3), with Weibull index £ = 2, we have 

P(cr th < o- ( ) + An cr f ) = 1 - exp [ - A o- f (2cr + A o- f )] (45) 

and obtain 

P nc (cro) = exp { - i- [1 - exp ( - A (2 + A ) cr 2 )] } (46) 

and 



^ = (1 - A () ) 



ln(l -A ) 



A„ (2 - A ) 



1/2 



(47) 



The corresponding results for arbitrary Weibull indices can be 
readily obtained. Pb(o~o) is again determined numerically by 
solving the recursion of Eqs. (37) and (38), with 



fi(o-„) = -^-{l - exp [ - A o-„ (2o- () + An cr„)]| . 



(48) 
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Figure 5: Ao-model with Weibull-distributed strength thresholds (k = 2). (a) 
No-cascade probability P nc and (b) breakdown probability Pb as a function of 
the initial stress ctq for different values of the load-redistribution parameter Ao. 
Lines: Results from Eqs. (46) [panel (a)] and Eqs. (37), (38), and (48) [panel 
(b)], respectively. Symbols: Results from Monte-Carlo simulations of the fail- 
ure process for N = 1000 fibers averaged over 1000 realizations. The statistical 
error is of the order of the size of the symbols. 



5.7.5. Discussion 

The behavior of P nc and 7\ as a function of crn and Ao is il- 
lustrated in Figs. 5 and 6 for a truncated Weibull distribution of 
strength thresholds, and we note that the results for uniformly 
distributed strength thresholds (see Sect. 5.1.1) show a quali- 
tatively similar behavior. Figure 7 shows the dependence of 
o-Q on Ao for both uniformly and Weibull distributed strength 
thresholds. 

In Fig. 5 and 6, our analytical results are compared with those 
obtained from Monte-Carlo simulations. To ensure the valid- 
ity of condition (2), we have chosen to keep Ao fixed as the 
number 7Vi n of intact fibers decreases, and to use the scaling 

po = Ao/7Vi n . 

To compare the results of this model (Ao-model) with those 
of the model analyzed in Section 4 (y-model), we first make the 
following observations. The y-model reproduces the GLS-limit 
if j < 2 and (in a particular sense) approaches an LLS-limit as 
y — > oo, while in the Ao-model, the GLS-limit is reproduced if 
Ao — > and the LLS-limit (transfer of the failing load to a single 
surviving fiber) for Ao — ♦ 1. Because of the different nature 
of the two A-distributions, an exact relation between y and Ao 
cannot be derived. A comparison of Figs. 2 and 7, however, 
suggests that for 2 < y < 6, a rough correspondence between 
the y- and the Ao-model is obtained if we set 



2/7 * 1 - Ao 



(49) 
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Figure 6: Ao-model with Weibull-distributed strength thresholds = 2). 
Breakdown probability as a function of the load-redistribution parameter Ao 
for different values of the initial stress ctq. Lines: Results from Eqs. (37), (38), 
and (48). Symbols: Results from Monte-Carlo simulations of the failure pro- 
cess for N = 1000 fibers averaged over 1000 realizations. The statistical error 
is of the order of the size of the symbols. 
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Figure 7: Ao-model. Critical stress cr^ as a function of the load-redistribution 
parameter Ao. Solid line: uniformly distributed strength thresholds [Eq. (43)]. 
Dashed line: Weibull-distributed (k = 2) strength thresholds [Eq. (47)] . 



Using this relation, it can be seen that the P nc vs. <x and P\, 
vs. cr behavior of the two models is qualitatively very similar 
if2<y<6(0<A < 0.7), except that in the y-model, P b (<T ) 
always increases towards one, while in the Ao-model, Pb(o"o) 
saturates at a value smaller than one if Ao > 0.1. 

For y > 6 (Ao > 0.7), however, the behavior of the two 
models is significantly different. The critical strength cr c Q , e.g., 
continues to decrease towards zero in the Ao-model, but starts 
to increase again in the y-model. Also the behavior of the no- 
cascade probability P nc vs. <t is completely different in the two 
models. In the y-model, P nc continues to increase towards one 
as y — > co, even for large values of <x , while in the Ao-model, 
P nc is bounded by P nc < 1 /e for large cr . As already discussed 
in Sect. 4.4.1, the peculiar behavior of the y-model for large y 
can be attributed to a breakdown of the continuum approxima- 
tion on which the corresponding A-distribution is based. 

6. Conclusions 

In this paper, we have introduced and analyzed a new fiber- 
bundle model with stochastic load redistribution. The fraction A 



of a failing fiber that is transferred to the surviving fibers is as- 
sumed to be a random variable, and we have considered two dif- 
ferent distributions for the A-values. The first (y-model) refers 
to a stochastic version of the range-dependent load redistribu- 
tion model of Hidalgo et al. [12], and the second (Ao-model) to 
a model with a simple bimodal A-distribution that can also in- 
terpolate between the two limiting cases of global and local load 
sharing. For the distribution of strength thresholds, we have 
also considered two different cases, a uniform and a Weibull 
distribution, both truncated below some finite stress o-q. 

While our models neglect any spatial correlations in the load 
redistribution after a failure, they have the advantage that they 
can be treated analytically, in contrast to most of the existing 
fiber bundle models that can only be analyzed via Monte-Carlo 
simulations. 

In the limit of global load sharing (y < 2 in the y-model 
or Ao — > in the Ao-model), our models not only recover the 
known exact results for the critical stress cr^, but also give the 
exact behavior of the breakdown probability for <t > cr^. In 
this GLS limit, the recursion relations for the determination of 
Pb(o"o) are reduced to simple transcendental equations 



Pbio-Q) = 1 - exp 



■ Pb(cr Q ) 



(50) 



1 - cr 

for a truncated uniform strength-threshold distribution, and 

P b (cro) = 1 - exp [ - 2al P b (cr Q )] (5 1 ) 

for a truncated Weibull distribution with index k — 2. Eqs. (50) 
and (51) can easily be derived from the recursion of Eq. (37) 
by taking the limit Ao — > 0, or from the transcendental equa- 
tion (28) for y < 2. 

With our stochastic models, we can also determine the crit- 
ical stress erg and the behavior of Pb(co) for <tq > cr^ in the 
case of a more localized stress redistribution (y > 2 in the y- 
model or Ao > in the Ao-model). As already discussed in 
Sect. 4.4.1, our y-model results for cr^ (for truncated uniform 
strength-threshold distributions) agree very well with the corre- 
sponding results of Raischel et al. [11] up to y w 6 or 7. It is 
quite remarkable that a stochastic model that neglects any spa- 
tial correlations can so accurately reproduce the behavior of a 
microscopically more adequate model. In addition, our analyt- 
ical solution allows us to trace back the onset of the transition 
between the GLS and LLS behavior at y = 2 to a change in the 
scaling of the upper cutoff A max of the A-distribution in the limit 
of infinite system sizes. 

In the case of strength-threshold distributions that are not 
truncated, the usual procedure is to gradually increase the ex- 
ternal force F from zero up to the complete breakdown of the 
entire bundle. The critical strength of the fiber bundle is then 
defined as the maximum stress F/N the system can support be- 
fore it breaks down. 

In global load sharing models, the surviving fibers always 
carry the same stress, so that the critical fiber-bundle strength 
can be written as 



°"bundle =Cr c[l -^th(0-c)]. 



(52) 



where cr c is the stress a surviving fiber carries at breakdown and 
1 ~ ^th(o"c) is the fraction of surviving fibers. This result is also 
recovered within our approach, with cr c = cr^. 

For nearest-neighbor LLS models, however, o-^ undle vanishes 
in the limit of large system sizes, o"£ undle oc 1 / y/N. As our 
stochastic models neglect spatial correlations, they cannot de- 
scribe such situations. 

We can, however, compare the results of our models with 
that of Ref. [15], where also a stochastic load redistribution 
model is used. Here, o"£ undle remains finite, even if the failing 
load is transferred only to a small, fixed number (n = 1,2,...) 
of randomly chosen surviving fibers. For n — 2 and for a 
uniform strength-threshold distribution, e.g., it is found that 
°"bundie ~ ^-2 f° r large systems. This can be compared with the 
corresponding result of our Ao-model. If we choose Ao = 0.5, 
so that the failing load is, on average, transferred to two surviv- 
ing fibers, we obtain cr^ = ( 1 - Ao)/(2 - Ao) = 1/3 for a uniform 
distribution of strength thresholds, i.e., 

<u„die = o%(l-o* ) = 2/9 * 0.22. (53) 

The discrepancy between the two models can be attributed to 
our assumption that at breakdown all surviving fibers carry the 
same stress cr^, whereas it is shown in Ref. [15] that at break- 
down, the surviving fibers have a broad distribution of stresses, 
with a pronounced exponential tail. 

It remains to be investigated whether our approach can be 
adapted to correctly analyze the behavior of fiber-bundle mod- 
els with a strength-threshold distribution that is not truncated. 
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